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A torsion creep theory for noncircular tubes was developed by applying a multiaxial creep theory 
to a derived generalization of Bredt's equations. 

A review of the literature revealed no torsion creep data on noncircular tubes. Hence, to evaluate 
the theory, a test program was carried out on twelve specimens of aluminum alloy structural tubing, 
of four configurations, at 400 °F. 

Observed discrepancies between the torsion creep theory and experiment are smaller than varia- 
tions in the measured creep properties of the specimen material from one tube configuration to another 
and are not appreciably greater than discrepancies between elastic torsion theory and experiment. 
Most of the observed discrepancies are consistent with measured anisotropy in the tubes, while other 
discrepancies are ascribed to nonhomogeneity in creep properties and a hydrostatic stress effect in 
multiaxial creep. 

For the calculation of torsion stresses in circular tubes the thin-wall approximation is adequate 
for thickness-to-radius ratios up to one-tenth. For straight-sided tubes equivalent accuracy is obtained 
for effective ratios up to only one-twentieth. These criteria apply to creep conditions as well as to elastic 
conditions. 
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1 . Introduction 

Torsion testing of thin-walled tubes furnishes a 
number of advantages over tensile testing in the de- 
termination of creep properties. The most significant 
of these is that the measurement of strain in torsion 
tests is essentially unaffected by the thermal expan- 
sions and contractions which accompany minor tem- 
perature fluctuations. Also, if the wall thickness is 
small compared to the diameter, the shear stress across 
the section may be regarded as uniform and, for a 
constant torque, this stress remains constant through- 
out the duration of the test. A final advantage, for cer- 
tain research purposes, is that it is comparatively easy 
to reverse the stress direction. 

In spite of these advantages most laboratory creep 
testing is carried out in tension. On the other hand, 
most structural members are subjected to multiaxial 
stress conditions. It is necessary, therefore, to relate 
creep behavior under multiaxial stress states to the 
uniaxial creep behavior that is usually obtained in the 
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laboratory. The purpose of this paper is to develop a 
means of calculating the torsion creep behavior of 
noncircular tubes from the tensile creep properties of 
the tube materials. 

The first recorded torsion creep tests of hollow tubes 
were carried out by Bailey [1, 2] 2 in 1929. Other 
pioneers in this area of research include Everett [3, 4] 
and Lea [5]. A. E. Johnson [6-10] has probably done 
more work in this field than anyone else and, indeed, 
with the exception of one related investigation [11], 
Johnson was the only researcher of record engaged in 
torsion creep studies between 1940 and 1960. Recently, 
however, there has developed a mounting international 
interest in the subject as evidenced by the appearance 
of six papers [12-17] in the last three years. This is 
attributable, in part at least, to anticipated creep 
deformation problems in the wings of supersonic 
aircraft. 

The tubular specimens used in the above investi- 
gations were, in every case, circular. This investi- 
gation is believed to be the first study of torsion creep 
in noncircular tubes. 

Most of the torsion creep tests cited in the literature 



2 Figures in brackets indicate the literature references on page 220. 
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were performed on miniature test specimens having 
diameters of about V2 in and gage lengths of about 2 
in. In fact, until recently [17] no torsion creep data 
had been reported on any tubular specimen longer 
than 8 in or larger than 1.6 in in diameter. 



2. Theory 

The theory of elastic torsion in thin-walled tubes, 
as developed by Bredt [18] in 1896 represents an 
approximation to Saint-Venant's general elastic torsion 
theory. Bredt's analysis was based on certain assump- 
tions regarding the nature of the stress distribution 
and results in two familiar equations. The first of these, 
which gives the stress distribution, is 



T 
2Ah 



(1) 



where 



r= shear stress 

T= torque 

A = area enclosed by the midthickness curve C of the 

tube cross section 
h = wall thickness. 

This equation is derived strictly from equilibrium con- 
siderations and is, therefore, assumed to be applicable 
to creep conditions as well as to elastic conditions. 

Bredt's second equation, which gives the angle of 
twist of the tube, is applicable only to elastic materials. 
However, a generalized version of Bredt's second 
equation is derived in appendix A and is applicable to 
material obeying any constitutive equation. This gen- 
eralized version of Bredt's second equation is 



= 



y m P 
2A 



(2) 



where 

6= angle of twist per unit length 
y m =mean shear strain 
P= perimeter of C. 

It is only necessary to insert the proper expression for 
y m in this equation to make it applicable to the par- 
ticular material behavior under consideration. 

For the present case, let it be assumed that the 
tensile creep behavior of the tube material can be 
expressed in the common form 



e = ko- m t" 



(3) 



where 



extensile creep strain 
ct= uniaxial tensile stress 
t = time 



and A, m, and n are constants which depend only on 
the material and the temperature. This equation 
applies only to constant stress, constant temperature 
conditions. 

If the octahedral shear stress theory of multiaxial 
creep is invoked, it can be shown that for a material 
obeying eq (3) the shear creep characteristics are 
given by 



y = V3A:(V3T) m * n 



(4) 



where 

7= shear creep strain 

and k, m, and n are the tensile creep constants defined 

by eq (3). Substituting eq (1) into eq (4) gives 



y 



*<$£)"*■ 



For the case of a thin-walled tube in torsion, a mean 
value for the shear creep strain may be obtained by 
performing a contour integration of this expression 
around C. That is, 



ym z 



\f-^§¥"\ 



ds 
h m 



(5) 



where 

g?s — incremental arc length of C. 

Equation (5) can now be inserted into eq (2) to give the 
angle of twist, 



2 A \ 2 A) Je h"< 



(6) 



This represents the creep counterpart of Bredt's sec- 
ond equation. If the tube in question has a uniform 
wall thickness h, eq (6) simplifies to 



, V3 P . VS T \ m 

1 = — -- k t n . 

2 A \ 2 Ah 



(7) 



For a circular tube with mean radius r m this equation 
may be written as 



r m \ 2 Ah) 



t" 



(8) 



Whereas eq (7) is approximate in the sense that it was 
derived with Bredt's assumptions regarding the nature 
of the stress distribution, eq (8) is exact in the sense 
that it coincides with a solution (appendix B) which 
does not require these assumptions. 
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3. Test Program 

A test program was designed to evaluate the appli- 
cability of eq (7) to several shapes of thin-walled tubes. 
Tension creep tests of the tube materials were carried 
out to determine k, m, and n, and torsion creep tests 
of the tubes were performed to provide experimental 
data for comparison with eq (7). 

3.1. Specimens 

Drawn seamless 2024-T3 aluminum alloy tubing 
was acquired for use as test specimens. This alloy 
was selected because of its reported adherence to 
the octahedral shear stress theory [19, 20] and because 
of its applicability to airframes. Four tubes of each of 
four different shapes were obtained: circular, square 
with rounded corners in two thicknesses, and rectangu- 
lar with rounded corners. The diameter of the circu- 
lar tubes, the sides of the square tubes, and the 
longer side of the rectangular tubes were all nomi- 
nally 3 in, and the lengths of the tubes were each 7 ft. 
The nominal dimensions of the four shapes are shown 
in figure 1. For any single tube, variations in the outside 
dimensions and in the thicknesses did not exceed 0.01 
in and 0.004 in, respectively. 
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FIGURE 1. Nominal dimensions of tube specimens, in inches. 

3.2. Tension Creep Tests 



14 data points from each test, a least squares approach 
was set up to find appropriate values of A:, m, and n for 
the material of each tube shape. The resulting ten- 
sion creep equations, in the form of eq (3), are as 
follows: For the circular tube material 



i3.4o- , -°V- :w X 10- 6 . 



For the 0.083-in wall square tube material 



(9a) 



e = 2.54o- , - 6 V , - 2H X 10- 6 . (%) 

For the 0.109-in wall square tube material 

e = 6.49o- 1 - 3 V- 36 Xl0- 6 . (9c) 



For the rectangular tube material 



e = 2.43cr 1 - 57 ^ 30 Xl0- 6 . 



(9d) 



In these equations cr is in kips/in 2 and t is in hours. 

A comparison between eq (9b) and the experi- 
mental tension creep curves for the specimens cut 
from the 0.083-in wall square tube is given in figure 
2. The discrepancies appear to be random and are 
attributed to the usual scatter in creep behavior. 
The discrepancies obtained with the specimens cut 
from the other tube shapes are similar to that shown 
in figure 2. 
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Figure 2. Tension creep behavior of specimens cut from 0.083-in 
wall square tube. 



Four tension creep specimens were cut longitudi- 
nally from one tube of each of the four shapes. Using 
a conventional dead weight, lever arm testing ma- 
chine, creep tests of approximately 50-hr duration 
were carried out on these specimens at 400 °F. The 
stress levels ranged from 9,100 to 22,500 lb/in 2 . Using 



To compare the tensile creep properties of the four 
tube materials, creep curves were calculated from 
eqs (9) for a stress of 15,000 lb/in 2 . These are given in 
figure 3. It may be seen that significant differences 
exist between the creep curves for the materials of 
two of the tube shapes as compared to those for the 
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Table 1. Torsion creep tests 
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FIGURE 3. Tension creep curves for the tube materials at 15,000 
lb I in 2 according to eqs (9). 

materials of the other two tube shapes. These dif- 
ferences may be attributable to variations in chemical 
composition or heat treatment, or to differences in 
the degree of .cold work resulting from the respective 
drawing operations. 

3.3. Torsion Creep Tests 

Twelve torsion creep tests were carried out at 
400 °F, three on each of the four tube shapes shown 
in figure 1. Each tube specimen was seven feet long 
and was mounted such that it passed symmetrically 
through a six-foot test furnace. Twist measurements 
were made over a three-foot gage length at the center. 
The twist measuring system had errors not exceeding 
one percent and individual measurements were re- 
peatable to better than 3 X 10~ 6 radians per inch. 
Temperatures were constant within 1 °F throughout 
each test and uniform, over the gage length, within 
3 °F. The testing equipment, the test procedure, and 
the calibration technique are described in appendix C. 

Table 1 lists the specimen designations, the torques, 
and the stress levels computed from eq (1) for each 
test. The experimental creep data are plotted in figures 
4, 5, 6, and 7 in terms of angle of twist against time. 
Also given in the figures are the theoretical torsion 
creep curves calculated with eq (7). For the circular 
tubes, figure 4, and the rectangular tubes, figure 7, 
agreement between theory and experiment is good for 
the lowest stress level. At the higher stress levels the 
theory predicts considerably less creep than was ob- 
served in the tests. For the 0.083-in wall square tubes, 
figure 5, the theory predicts less creep at all stress 
levels than was actually obtained, while for the 0.109- 
in wall square tubes, figure 6, the theory predicts 
significantly more creep at all stress levels than was 
actually obtained. Although the relationship between 
experiment and theory is reasonably consistent for 



Specimen 

No. 


Torque 


Shear 
stress 


Stress 
error 




in-lb 


lb/in 2 


Percent 



Circular tubes 



1 


9,930 


8,940 


+ 2 


2 


11,800 


10,600 


-25 


3 


10,900 


9,820 


-23 



0.083-in wall square tubes 



4 


11,800 


8,080 


-15 


5 


13,800 


9,380 


-18 


6 


15,600 


10,700 


-12 



0.109-in wall square tubes 



7 


16,600 


8,920 


+ 31 


8 


18,500 


9,940 


+ 30 
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14,700 


7,900 


+ 23 



Rectangular tubes 
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FIGURE 4. Torsion creep behavior of circular tubes. 

each tube shape, no dominant relationship is apparent 
from all of the data taken together. 

The discrepancies between calculated and experi- 
mental results were evaluated in terms of stress. That 
is, for each test the stress that would have brought 
the calculated curve into approximate coincidence 
with the experimental data was computed. The dif- 
ference between this value and the theoretical stress, 
expressed as a percentage, is given in table 1 as the 
"stress error." 
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FIGURE 5. Torsion creep behavior o/0.083-in wall square tubes 
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FIGURE 6. Torsion creep behavior of 0.109-in wall square tubes. 
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Figure 7. Torsion creep behavior of rectangular tubes 
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4. Analysis of Results 

The discrepancies between the theoretical and 
experimental creep curves are not significantly 
greater than the discrepancies that are obtained 
when Bredt's equations are applied to elastic tubes at 
room temperature [21-24]. This, however, does not 
explain the reasons for the discrepancies and, there- 
fore, gives no indication of the limits of applicabil- 
ity of the theory. 

In order to explain the observed discrepancies be- 
tween theory and experiment, it is desirable to re- 
examine the justifications for the numerous assump- 
tions which were involved in the theory. These assump- 
tions may be divided into three categories. The first 
category includes those common assumptions regard- 
ing the idealization of the specimen materials which 
are necessary to make them amenable to mathematical 
treatment. These assumptions are: 

1. The material is continuous. 3 

2. The material is homogeneous. (See footnote 3.) 

3. The material is isotropic. 

In the next category are the additional assump- 
tions inherent in the Saint-Venant torsion theory. 

4. Ends are free to extend or contract in the length- 
wise direction. 

5. Ends are loaded by a pure torsional couple with 
no resultant force. 

6. There is no restraint of warping. 

7. There are no reentrant corners or other areas of 
severe stress concentration. 

Finally, there are the assumptions required by Bredt's 
approximation to the Saint-Venant theory. 

8. The shear stress parallel to the surface is uniform 
throughout the wall thickness. 

9. The shear stress normal to the surface is zero 
throughout the wall thickness. 

4.1. Material Idealization 

a. Continuity 

Although it is known that metallic materials are 
composed of crystals, their dimensions are generally 
very small in comparison with the dimensions of the 
specimens. Hence, the mathematical treatment of the 
material as though it were continuous is generally 
reasonable. In the present case, however, there was 
some doubt as to whether the crystals were, in fact, 
small in comparison with the wall thicknesses of the 
tubes. To resolve this doubt, small samples were cut 
from one square tube of each of the two wall thick- 
nesses. For each sample the number of grains in the 
thickness direction was counted both on a longitudinal 
surface and on a transverse surface. In all cases the 
count ranged between 33 and 48 grains. On this basis, 
it is felt that the assumption of continuity is justified. 



3 For these two characteristics, continuity and homogeneity, some authors prefer the 
terms homogeneity and uniformity, respectively. In the present context continuity indi- 
cates the nature of a structureless mass and homogeneity is the quality of having the 
same mechanical properties throughout the material, i.e., independent of location. 
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b. Homogeneity 

It is conceivable that the tubes from which the ten- 
sile creep specimens were cut may not have been 
representative of the tubes which were tested in tor- 
sion creep. Therefore, the homogeneity of the specimen 
materials was checked, as follows, with a hardness 
survey. 

A rectangular plate, approximately 1 in by 2 in, 
was cut from each torsion creep specimen, from the 
region which had been between the grips and the 
test furnace. The material in this region is consid- 
ered to have suffered the least change in proper- 
ties due to stressing or heating. In addition, a similar 
plate was cut from each of the four tubes from which 
the tension creep specimens had previously been 
cut. 

Table 2 lists the average Vickers hardness numbers 



Table 2. Hardness test i 



<ilts 



Specimen 


Hardness 




No. 








VHN 




1 


146 




2 


148 




3 


141 




C 


143 


Average. 




144 




4 


147 




5 


148 




6 


145 




SN 


147 


Average. 




147 


7 


145 




8 


148 




9 


142 




SK 


146 


Average. 




145 




10 


145 




11 


140 




12 


141 




R 


140 


Average. 
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for each of the test plates. The plate numbers in the 
first column refer to the tube specimen numbers 
(table 1) from which the respective plates were cut. 
The plates designated by letter symbols refer to 
those tubes from which the tensile creep specimens 
had been cut. It can be seen that the variation in 
hardness within each batch is small, as expected. 
Unfortunately, the data also indicate that the dif- 
ference in hardness from batch to batch is also small. 
Yet, it has been determined (fig. 3) that a signifi- 
cant difference in creep properties exists between 
different batches. Hence, it must be concluded that 
the uniformity in hardness within each batch does 
not necessarily indicate a uniformity in creep prop- 
erties. The possibility thus remains that the dis- 
crepancies between the calculated and experimental 
torsion creep behaviors may be due to nonhomoge- 
neity. 



c. Isotropy 

As a result of the cold-working process, drawn tubes 
are typically orthotropic with the maximum strength 
properties coinciding with the longitudinal direction. 
In the present investigation the tension creep prop- 
erties of the tube materials were determined from 
specimens cut longitudinally from the tubes. In the 
torsion tests, however, the principal directions are at 
45 deg to the longitudinal axis of the tubes. It is nat- 
ural, therefore, to question whether the creep prop- 
erties of the material in the diagonal (45-deg) direction 
and in the longitudinal direction are equivalent. Un- 
fortunately, the tensile creep properties in the diagonal 
direction could not be evaluated in the present investi- 
gation because the specimen lengths which could be 
cut diagonally are too short to be satisfactorily ac- 
commodated by the available creep testing equipment. 

As a compromise, it was decided to secure a general 
picture of the anisotropy in the tubes by conducting 
room-temperature tensile tests on miniature speci- 
mens which had been cut from a tube in the longi- 
tudinal, diagonal and transverse directions. The first 
set of miniature specimens was taken from the 0.083- 
in wall square tube from which the tension creep 
specimens had been cut. In figure 8 the longitudinal, 
diagonal and transverse stress-strain characteristics 
are compared. It may be seen that anisotropy does, 
indeed, exist with the material being stronger in the 
longitudinal direction than in the transverse direction. 

Tensile and yield strengths from these tests are 
given in table 3, each tabulated value being the aver- 
age from a pair of tests. The ratio of yield strengths 
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Figure 8. Comparison of tensile stress-strain characteristics of 
0.083-in. wall square tube material in three directions. 
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TABLE 3. Tensile properties of 0.083-in wall square tube material I.6XI0" 3 



Direction 


Yield 

strength ;| 


Tensile 
strength 


Longitudinal 


Iblui- 
52,000 
45,200 
43,500 


Iblin 2 
72,800 


Diagonal 


69,200 


Transverse 


67,800 







a 0.2 percent offset. 

in the longitudinal and diagonal directions is 1.15. 
If this ratio may be taken as a measure of the plastic 
anisotropy, then an estimate can be made of the effect 
of anisotropy on the torsion creep behavior of the« 
0.083-in wall square tubes. On the basis of this ratio 
it is postulated, for purposes of this estimate, that 
plastic flow in the diagonal direction may be deduced 
from the plastic flow properties in the longitudinal 
direction by considering the stress in the diagonal 
direction to be 15 percent greater than the true value. 
Thus, the creep strain in the diagonal direction 4 is 
taken as [cf eq (3)| 

€ = /-(1.15o-)"7". 

With this modification the torsion creep equation, 
eq (7), becomes 



V3 P , /V3 1.157Y" 
6 = k\ t". 

2 A V 2 Ah 



(10) 



An independent justification for this form of aniso- 
tropy correction is given in appendix D. Torsion creep 
curves were calculated for the 0.083-in wall square 
tubes using this equation. These calculated curves 
are compared with the experimental results in figure 9. 
It is seen that, on the whole, the agreement between 
calculated and experimental results is now quite 
acceptable. The remaining discrepancies appear to 
be random, rather than systematic, and are apparently 
within the experimental variability of the creep prop- 
erties of the material as exhibited by figure 2. 

On similar bases it is reasonable to presume that the 
discrepancies obtained with the circular and rectangu- 
lar tubes (figs. 4 and 7) can likewise be explained, for 
the most part, in terms of anisotropy. Unfortunately, 
however, this type of correction would increase the 
discrepancies obtained with the 0.109-in wall square 
tubes (fig. 6) since, in this case, the isotropic theory 
predicts considerably more creep than was observed 
in the tests. This situation naturally creates a sus- 
picion that the 0.109-in wall square tubes possess a 
reverse type of anistropy: that is, one where the mate- 
rial is stronger in the diagonal (and transverse) direc- 
tion than in the longitudinal direction. 
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Figure 9. Torsion creep behavior of 0.083-in wall square tubes. 

Comparison of experimental data with theory corrected for anisotropy. 

To check this, a set of miniature specimens was cut 
from the 0.109-in wall square tube from which the 
tension creep specimens had earlier been taken. The 
stress-strain characteristics determined from these 
specimens indicate that the anisotropy in this case is 
qualitatively similar to that obtained with the thinner 
tube material (fig. 8). It must be concluded therefore, 
that anisotropy cannot, by itself, explain all of the 
discrepancies between experimental and calculated 
torsion creep behavior, and that some other factor or 
factors, as yet unknown, are operative here. 

The tensile and yield strengths for the 0.109-in wall 
square tube material are given in table 4. 

Table 4. Tensile properties of 0.109-in wall square tube material 



Direction 


Yield 
strength a 


Tensile 
strength 


Longitudinal 


lb/in 2 
55,700 
47,200 
43,300 


lb/in 2 
70,900 
69,600 
67,500 


Diagonal 


Transverse 





4 The degree of anisotropy is so small that the angle between the directions of principal 
stress and principal strain may be ignored. 



a 0.2 percent offset. 

d. Compressive Creep Properties 

Another assumption of concern here, which is 
related to the matter of isotropy, is the equivalence of 
the creep properties in tension and in compression. 
This was not listed as a separate assumption but it 
is inherent in the octahedral shear stress theory of 
multiaxial creep, which was used in the development 
of the torsion creep equations. Vawter et al. [25], 
and Heimerl and Farquhar [26] measured and com- 
pared the tensile and compressive creep properties of 
2024-T3 aluminum alloy at several temperatures. 
Although no data were obtained at 400° F their re- 
sults suggest that at this temperature the creep prop- 
erties in tension and compression are probably close 
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to one another. Similar results have been reported 
[27] for 2024-T4 aluminum alloy. 

4.2. End Effects 

Assumptions 4, 5, and 6 may be considered collec- 
tively under the heading of end effects. Conformance 
with these three assumptions is required in order to 
achieve, in the specimen, the state of pure torsional 
loading which is considered in the theory. 

a. Longitudinal Restraint 

If a tubular specimen were to manifest a tendency 
to lengthen or shorten under torsion, and if this tend- 
ency were resisted, axial stresses would be induced 
in the tube. Accordingly, the design of the torsion 
creep testing equipment incorporated provision for 
the torquing shaft to slide axially in its bearings. 
Thus, with exception of some friction, the torqued end 
of the specimen was free to extend or contract. 

b. Applied Stress Distribution 

The Saint-Venant solution to the torsion problem 
concludes that the stress distribution is the same at 
every cross section. This indicates that, for conform- 
ance with the theory, the torque must be applied to 
the end surface of the cylinder as a system of trac- 
tions identical with the shear stress distribution which 
is thence obtained at every cross section. The restric- 
tive nature of this requirement is eased by the appli- 
cation of Saint-Venant's principle, which is usually 
satisfactory. In the case of thin-walled members, 
however, it has been known for some time that Saint- 
Venant's principle is not universally valid. 

With a thin-walled tube in torsion, if the applied 
system of shear forces is not identical with the derived 
shear stress distribution, then longitudinal normal 
stresses are created in order to maintain compatability 
of the deformation at the ends. (This is illustrated 
later for the case of a square tube.) These normal 
stresses differ from those induced by the restraint 
of elongation or contraction in that they are self- 
equilibrating at any cross section. Nevertheless, these 
stresses may make themselves felt at longitudinal 
distances which are many times the depth or diameter 
of the tube. 

Thus, assumption No. 5, which is adequate for solid 
cylinders is inadequate for thin-walled tubes. Instead, 
strict conformance with the Saint-Venant torsion 
theory requires that the torque be applied in such 
a way that the applied distribution of shear forces is 
identical with the derived distribution of shear stresses. 
In a practical sense, deviations can be tolerated only 
if the induced longitudinal stresses are comparatively 
small. 

The dissipation of the longitudinal normal stresses 
with distance from the ends was calculated by Vlasov's 
method [28] for the square and rectangular tubes of 
the present investigation. The results of the calcula- 
tions are given in figure 10, which shows the variation 
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FIGURE 10. Dissipation of longitudinal normal stresses. 

of the longitudinal normal stresses a as a fraction of 
the indeterminate value cr at the ends. It may be seen 
that, owing to the fact that 21-in lengths of tube were 
provided between the grips and the gage length, the 
longitudinal stresses within the gage length have dis- 
sipated to less than five percent of the values at the 
ends of the tubes. This indicates that the conversion 
of the applied shear stress distribution into the Saint- 
Venant shear stress distribution is essentially ac- 
complished. Only in an infinitely long tube would this 
conversion be complete. 

The absence of significant longitudinal normal 
stresses within the gage length was checked experi- 
mentally with strain gages, as described in appendix C. 

c. Warping 

In the case of rectangular tubes, longitudinal nor- 
mal stresses are also introduced by the restraint of 
warping. The effects of these stresses on torsional 
stiffness were calculated by von Karman and Chien 
[29]. From their results it was determined that, for 
the rectangular tubes tested in the present investi- 
gation, the restrained warping causes a noticeable 
effect on the angle of twist for distances of only 
an inch or two from the grips. Thus, within the gage 
length, the warping may be considered to be unre- 
strained, as required by assumption No. 6. 

4.3. Stress Concentration 
The square and rectangular tubes used in the pres- 
ent investigation were selected with rounded corners 
in order to minimize the stress concentrations. These 
stress concentrations were evaluated from Huth's 
charts [30] as follows: 

For the 0.083-in wall square tubes, k c — 1.02; 
for the 0.109-in wall square tubes, k c = 1.15; 
for the rectangular tubes, k c = 1.00. 
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It may be seen that the stress concentrations are, 
indeed, small. The largest value, for the 0.109-in 
wall square tubes, cannot reasonably explain the dis- 
crepancies observed between the theoretical and ex- 
perimental torsion creep results. If anything, the stress 
concentration would tend to cause greater experi- 
mental deformations than the theory predicts whereas, 
in fact, the opposite was observed with the 0.109-in 
wall square tubes. 

4.4. Thickness Effects 

Assumptions 8 and 9, which are required by Bredt's 
analysis, may be discussed under the subject of 
thickness effects since, to a large extent, these 
assumptions are justified if the wall thickness is small 
enough. 

a. Stress Gradient 

Curiously, there is no well-known criterion or rule- 
of-thumb which defines the maximum thickness or 
other limiting condition for application of thin-wall 
torsion theory. On the other hand, for pressure ves- 
sels there is a widely quoted criterion that defines 
a thin cylinder as one in which the wall thickness is 
less than one-tenth of the radius. It is interesting to 
examine this criterion with a view toward establishing 
a similar one for torsion. 

(1) Circular Tabes. In an internally pressurized cir- 
cular tube the tangential stresses vary through the 
wall thickness, reaching a maximum at the inside 
surface. Under elastic conditions the maximum stress, 
as given by Lame's equation, may be written as 



-ax /'L+4 



where p is the pressure and a is the ratio of the wall 
thickness h to the mean radius r m . For thin-walled 
tubes the tangential stress is considered to be uniform 
through the thickness with a value of 



pn (I 1\ 



where n is the inside radius. The ratio 



imum at the outside surface. The maximum stress, 
as given by Saint-Venant's theory, may be written as 
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where T is the torque. For thin tubes the shear stress 
is considered to be uniform through the thickness 
with a value 
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The ratio of the maximum shear stress to the shear 
stress given by the thin-wall theory is 



l+o 

7"max ^ 



' max 
T thin 



1+ ? 

4 



(12) 



Equations (11) and (12) are plotted in figure 11. It is 
seen that the error involved in the thin-wall approxi- 
mations is greater for the pressure case than for the 
torsion case. Hence, a criterion applicable to the 
pressure case is reasonable for the torsion case as 
well. When a = 0.1 the error is slightly more than five 
percent for the pressure case and slightly less than 
five percent for the torsion case. For the circular tubes 
tested* in the present investigation, a = 0.057, which 
corresponds to an error of less than three percent. 
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is an index of the error of the thin-wall approximation 
for any specific value of a. 

For elastic torsion of a circular tube the shear 
stresses vary through the thickness, reaching a max- 



FlGURE 11. Effect of wall thickness on the thin-wall approximation 
of the stress. 

(2) Straight-Sided Tubes. For elastic torsion of 
square or rectangular tubes the situation is not as 
favorable. If a mean radius for straight-sided tubes is 
defined by 

2A 
r m = — (13) 

and if the stress concentrations at the corners are 
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neglected, then the maximum shear stress at the 
outer surface is [24, 31] 



2Ah 



fl + a) 



Therefore, 



Tthii 



-=l + a. 



Using the criterion of five percent error that was 
used above it appears that, for straight-sided tubes 
in torsion, thin-wall theory is applicable only if a is 
less than 0.05, rather than 0.1, which is the case for 
circular tubes. The following values of a apply to 
the square and rectangular tubes tested in the present 
investigation: 

For the 0.083-in wall square tubes, a = 0.057; 
for the 0.109-in wall square tubes, a = 0.074; 
for the rectangular tubes, a = 0.066. 

Clearly, the wall thicknesses for these tubes exceed 
the limit within which thin-wall torsion theory is ac- 
curate. However, this limit applies only to the stresses. 
It has been shown [32, 33] that the thin-wall require- 
ment is far less restrictive for the calculation of 
angle of twist than it is for the calculation of stresses. 
Hence, the use of the thin-wall approach is acceptable 
for present purposes. 

(3) Effect of Creep. For the case of circular tubes 
the exact creep stress distribution is given by eq (B.3) 
in appendix B. As an illustration of this stress dis- 
tribution, this equation, with r m =l and A=0.1, 
is plotted in figure 12 for creep exponents of m=l, 
m=2, and m=^. The curve for ra=l represents 
the elastic stress distribution. The curve for m — 2 
represents a creep stress distribution where elastic 
strains are negligible in comparison with the creep 
strains. For a real case where a combination of elas- 
tic and creep strains exist the distribution would He 
somewhere between the curves for ra=l and ra = 2. 
The curve for ra = °c represents the uniform stress 
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Figure 12. Variation of shear stress through the wall of a circular 
tube in torsion. 



distribution assumed by Bredt's thin-wall theory. 
Thus, it is seen that in creep the stress distribution is 
actually closer to the thin-wall approximation than it 
is under elastic conditions. It may be concluded, 
therefore, that any criterion which defines a limiting 
wall thickness for elastic torsion is acceptable for 
creep conditions, too. 

b. Stress Direction 

Assumption No. 9 requires that all shear stresses be 
parallel to the surface of the tube. Thinking in terms 
of the hydrodynamic analogy, this appears to be a 
reasonable assumption for tubes of smooth contour, 
such as circles, but it can be demonstrated that it 
is somewhat erroneous for tubes with corners. 

Consider, for example, torsion of a thin-walled tube 
of length L having a cross section in the form of a 
square with sides of length it and thickness h. This 
is shown diagramatically in two views in figure 13(a). 
According to the assumption, the shear stresses r in 
each side of the square are parallel to the surface. If 
these are considered to be the only stresses acting, 
the only deformation which could take place is pure 
shear with displacements 8 as shown in figure 13(b). 
(Only one face is shown in the left view, for sim- 
plicity. The broken lines indicate the original, unde- 
formed shape.) The angle y is simply the shear strain 
rjG\ where G is the modulus of rigidity and ihe prime 
is added for generality; G' and y may be considered 
time-dependent to allow for creep. Displacement 8 
is equal to yL. 
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FIGURE 13. Deformation of a square tube in torsion. 
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This deformation pattern is clearly unacceptable 
by itself. What is required to restore compatability to 
the tube shape is a torsion of each face about its own 
longitudinal axis, as indicated by the torques T\ in 
figure 13(c). The final configuration is shown in figure 
13(d), where the required angle of twist is 
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The torques T\ are comprised of shear stresses near 
the corners which act normal to the surface although 
they are necessarily zero at the surface. The magni- 
tude of T\ is simply the torque required to twist a flat 
strip through angle i//. That is, 



7\ = 



G u 3 



r/r 3 



where J is the torsional constant. 

The condition of moment equilibrium can now be 
imposed. The four torques Ti plus the moment of the 
shear stresses r [fig. 13(a)) must equal the applied 
torque T. That is, 

T=W x + ^Thu)- = -Th* + 27lui 2 

from which 



2h l-h 2 + u 2 



According to Bredt's equation, however, the shear 
stress is 
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The difference, expressed as a fraction of the Bredt 
stress is 
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This quantity is negligible for the square tubes tested 
in the present investigation and, in fact, amounts to 
less than one percent for hlu ratios up to one-tenth. 
The presence of shear stresses normal to the sur- 
face of the tube introduces another source of error, 
namely, distortion of the cross-sectional shape in its 
own plane. This distortion changes the enclosed 
area and, therefore, the stresses as well. This prob- 
lem was studied in detail by Ikeda [34] for square, 
thin-walled tubes. Using the results of his analysis 
it was determined that the effect on the stresses is 
less than one percent for the square tubes tested in 
the present investigation. 



5. Discussion 

The foregoing analyses have demonstrated that 
the assumptions introduced by Bredt's approxima- 
tion to the Saint-Venant torsion theory are appli- 
cable to the torsion creep problem studied in this 
investigation. In the area of material idealization, 
however, some aspects of questionable applicability 
were revealed. Specifically, it was found that the 
materials of the tubes tested were not isotropic. 
However, with an approximate correction for anisot- 
ropy, satisfactory coincidence between theory and 
experiment was shown to be attainable. Unfortu- 
nately, this was not the case for the 0.109-in wall 
square tubes. 

The studies carried out were inadequate to demon- 
strate that the specimen materials are homogeneous 
insofar as creep properties are concerned. There- 
fore, the discrepancies between theory and ex- 
periment could logically — although tentatively — be 
ascribed to nonhomogeneity; little can be done with 
such a problem unless the actual distribution of 
creep properties in the specimen is known. Intuitively, 
however, this is an unsatisfactory conclusion. The 
discrepancies (fig. 6) appear too large and too con- 
sistent to be attributed to a negative characteristic 
such as an absence of homogeneity. The fact that the 
discrepancies are not appreciably greater than those 
obtained under elastic conditions does not mitigate 
the need for a more satisfactory explanation, and at 
least one other possible explanation does exist. 

In the preceding sections most of the assumptions 
involved in the development of the theory were justi- 
fied or evaluated by measurements or calculations 
applicable to the problem at hand. The justifications 
for two of the assumptions, however, were based on 
references to data in the literature. These two as- 
sumptions, which are related, are the octahedral 
shear stress theory of multiaxial creep, and the 
equivalence of tensile and compressive creep. The 
failure of the torsion creep theory to describe the be- 
havior of the 0.109-in wall square tubes is reason 
enough to reexamine the adequacy of these justi- 
fications. Although the octahedral shear stress theory 
is reported to be the one most commonly used at the 
present time [35] and also most favored by the weight 
of evidence in the literature [36J, careful experimental 
measurements have been presented, in recent years, 
which are inconsistent with it [12, 37]. 

Some investigators are inclined to the belief that 
inconsistencies with the octahedral shear stress theory 
are somehow attributable to an effect of hydrostatic 
stress. The theory neglects such an effect but a dif- 
ference between tensile and compressive creep prop- 
erties would indicate that the effect is not negligible. 
Finnie [12 1 attempted to evaluate this effect experi- 
mentally and found that it could not be characterized 
by any simple modification of the octahedral shear 
stress theory. By means of a literature search, he also 
concluded that the hydrostatic stress effect becomes 
significant only when the absolute test temperature 
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exceeds half the absolute melting temperature of the 
test material, and proposes this as an explanation for 
the belief (which was widespread until recently) that 
the hydrostatic stress effect is negligible. In this con- 
text it is significant that the tests in the present in- 
vestigation were conducted at a temperature which 
does, in fact, exceed half the absolute melting tempera- 
ture of the specimen material. 

The current state of knowledge regarding creep 
under multiaxial stresses was adequately summed 
up by Manson [35] in his discussion of the octahedral 
shear stress formulation: 

Whether the quantitative discrepancies that 
occur . . . are due to the inadequacy of the 
formulation still remains to be established. If 
modifications are to be made, the question must 
still be resolved as to whether they should be 
in the form proposed by Wahl [38], or by a suit- 
able correction which includes hydrostatic stress, 
or whether still another interpretation is best. 
It appears that an improved torsion creep theory for 
noncircular tubes must await the development of an 
improved multiaxial creep theory. 

6. Conclusions 

A theory for the torsion creep behavior of thin-walled 
tubes was developed by applying the octahedral shear 
stress theory of multiaxial creep to a generalized ver- 
sion of Bredt's equations. The theory was compared 
with experimental results obtained by tests on 12 
large circular and noncircular aluminum alloy tubes 
at 400 °F. The following principal conclusions, appli- 
cable to the material and tube configurations tested 
herein, were obtained: 

1. In most cases the theory can provide usable 
results when suitable allowances or corrections are 
made for anisotropy. 

2. Aside from the effects of anisotropy, differences 
between theory and experiment have been tentatively 
ascribed to (a) a shortcoming in existing theory of 
multiaxial creep behavior which is related to the 
inequality of tensile and compressive creep properties, 
or (b) nonhomogeneity in creep properties. 

3. The largest differences between theory and ex- 
periment were obtained with the 0.109-in wall square 
tubes. These differences were conservative and are 
smaller than the differences in the measured creep 
properties of the material from batch to batch. Also, 
the differences between theory and experiment are 
not appreciably larger than the discrepancies obtained 
when Bredt's equations are used to predict the tor- 
sion behavior of elastic tubes at room temperature. 

Other conclusions of interest which were obtained in 
the course of the investigation are: 

1. The existence of nonhomogeneity in creep prop- 
erties at elevated temperatures is not, necessarily, 
revealed by measurements of other mechanical prop- 
erties at room temperature. 

2. For the calculation of stresses in circular tubes 
the thin-wall approximation may be considered ade- 



quate for thickness-to-radius ratios up to one-tenth. 
For straight-sided tubes equivalent accuracy is ob- 
tained for ratios up to only one-twentieth, where the 
effective radius is defined as in eq (13). These cri- 
teria are applicable under creep conditions as well as 
under elastic conditions. For the calculation of twist 
angles the thin-wall approximation presumably is 
less restrictive than it is for the calculation of stresses. 
3. A simple and unique analysis is presented for 
evaluating the shear stresses which act normal to 
the surface of square tubes under torsion. It is shown 
that the effect of these stresses, on the overall stress 
distribution, is negligible for thin-walled tubes. 
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Appendix A. Generalized Version of Bredt's 
Second Equation 

Consider a cylindrical shell of closed section having 
a length L and faces normal to the cylinder axis. See 
figure 14. Under an applied torque T an originally 
longitudinal element AB of the midsurface rotates to 
an orientation AB' . Angle BAB' is the shear strain y. 
Displacement BB' is equal to Ly and it is assumed that 
the strains are sufficiently small that this displacement 
is perpendicular to the radius r from the axis of 
torsion. 

The angle of twist, viewed from this axis, is 



BB' = Ly 

r r 



(A.l) 
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hteratury, an d t ne an gl e f twist per unit length is 
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It is clear that this quantity does not remain constant 
as midthickness curve C is traversed. Yet it is com- 
mon and desirable to assign a single value to the 
quantity known as angle of twist per unit length. Gen- 
erally, this question is ignored and is computed either 
by the membrane analogy or by a strain energy ap- 
proach, both of which yield an effective mean value 0. 
A mean value is also sought here, but the membrane 
analogy and the strain energy method will not be used 
since these approaches apply principally to elastic 
materials. 

A mean radius can be deduced from the relation 
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where P is the perimeter of C and A is its enclosed 

L- 




FloiJRE 14. Derivation of generalized version of Bredt's second 
equation. 
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area. Then, designating the mean shear strain as y m 
leads to the following expression for angle of twist. 



y m 
r m 



JmP 

2A ' 



(A.4) 



This equation may be regarded as a generalized ver- 
sion of Bredt's second equation. 

Appendix B. Torsion Creep of Circular Tube 

Consider a hollow circular rod in torsion with inside 
radius a and outside radius b. The shear creep prop- 
erties of the material are given by eq (4). Then, ne- 
glecting elastic strains for the moment, 5 the shear 
stress at any radius r is shown [39, 40] to be 
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The shear strain distribution is obtained by substi- 
tuting eq (B.l) into eq (4), giving 
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Then the angle of twist is 
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For a thin-walled circular tube let the mean radius 
be r m and the wall thickness be h so 
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Then the stress distribution, eq (B.l), becomes 
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Since - — <^ 1, the quantity in brackets in the denomi- 

/ 1 \ h 
nator reduces to approximately 2 ( 3 H — I- — . Thus, 

\ m) 2r m 

the stress distribution may be written as 

= T ( r v /m - 

27nj n h\r m ) 

5 The effect of the elastic component is discussed in section 4.4a(3). 
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Using a similar approach, the angle of twist, eq (B.2), 
may be rewritten as 
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This is identical with eq (8) which was derived with the 
generalized Bredt equation. Note also that the shear 
stress at the mean radius, from eq (B.4), is 



m 2nr l m h 2Ah 

which is the value of the uniform shear stress accord- 
ing to Bredt's first equation. 

Appendix C. Details of Torsion Creep Tests 
C.l. Equipment 

An overall view of the torsion creep testing equip- 
ment is shown in figure 15. Close-up views of the left 
and right ends of the setup are given in figures 16 and 
17, respectively. 

The ends of the specimen (A, figs. 15 and 16) are 
reinforced with close-fitting plugs to a depth of three 
inches. The specimen is mounted such that it passes 
through the test furnace (B, figs. 15, 16 and 17). The 
plugged left end of the specimen is rigidly bolted to 
a welded lug on a building column. The right end of the 
specimen is fastened, by means of a keyed shaft 
passing through the end plug, to a large sprocket 
wheel (C, figs. 15, 16 and 17). The shaft is supported 
by three journal bearings used at only a small fraction 
of their load capacity. The horizontal positions of the 
bearings are adjustable, both longitudinally and lat- 
erally, to permit proper alinement of the shaft with 
the fixed end mounting. The sprocket wheel, which 
has a pitch diameter of 38.2 in, applies torque to the 




FIGURE 15. Overall view of torsion creep testing setup. 
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View from IH 



test specimen by means of a weight pan (D, fig. 17) 
suspended from the sprocket chain. 

The test furnace is fabricated of transite, with 
steel reinforcement at the edges. The side of the 
furnace which faces the wall in figure 15 is hinged. 
Figure 18 is a view of the furnace, in a different 
location, with the hinged side opened to reveal the 
interior. Heat is supplied by 32 tubular infrared 
lamps rated at 120 V and 500 W. The lamps are wired 
in parallel pairs, the voltage to each pair being ad- 
justable by one of a bank of variable autotransformers 
(E, fig. 15). Power is regulated by three time-propor- 
tioning, potentiometer-type temperature controllers 
(F, fig. 15), each controlling the temperature over 
approximately one-third of the furnace. Nine chromel- 
alumel thermocouples were spot welded to the outer 
surface of each torsion specimen, at various locations 
within the three-foot gage length. Three of these 
thermocouples served as the sensors for the tempera- 
ture controllers and their outputs were continuously 
recorded. Outputs of the other six thermocouples 
were monitored with a calibrated millivolt potentiom- 
eter (G, fig. 15), which was also used to periodically 
check the calibration of the temperature recorders. 

Thermocouples cannot normally be spot welded 
to the structural aluminum alloys with portable, 
laboratory-type spot welding apparatus. However, 
the usual alternative mounting technique, clamping, 
was considered to be particularly undesirable for 
this type of specimen since the clamps would have to 
be large as well as numerous and, therefore, dele- 
terious to the attainment of a uniform temperature 
distribution. Hence, a new technique was devised 
for spot welding thermocouple beads to 2024-T3 



FIGURE 17. Torsion creep testing setup. 

View from right side. 



aluminum alloy, which turned out to be completely 
successful. The technique, in brief, consists of spot 
welding a small piece of 0.004-in gold foil, about 1 
mm square, to the aluminum and then spot welding 
the thermocouple bead to the foil. 

Twist was measured with a troptometer 6 developed 
specifically for this investigation. Steel sprockets 
(H, fig. 18) having a pitch diameter of 5.73 in were 
concentrically attached to the specimen at each end 
of the three-foot gage length. Sprocket chains (I, 
figs. 16 and 18) passed over these sprockets and out 
of the furnace through small holes and thence over a 
pair of small sprockets (J, figs. 16 and 17) having a 
pitch diameter of 0.81 in. Each small sprocket is 
mounted on the shaft of a precision, one-turn, resist- 
ance potentiometer rated at 10,000 fl ± 1 percent 
with a linearity tolerance of ± 0.50 percent. The ends 
of the chains are weighted (K, figs. 16 and 17) to keep 
them taut. 

Torsion of the specimen causes rotation of the 
5.73-in sprockets mounted on it. This rotation is 
transferred, by means of the sprocket chains, to the 
resistance potentiometers, and magnified by the 
sprocket ratio. The difference in rotation of the two 
5.73-in sprockets, measured as the difference in resist- 
ance between the two potentiometers, is the angle 
of rotation of the specimen over the three-foot gage 

* tn5fp.tom'"e".ter, n. An instrument for measuring the angular distortion of a bar or piece 
undergoing a torsion test (Webster). From the Greek tropos, turn and metron, measure. 
This word is clearly preferable to any of the existing words used for this purpose in the 
current literature. 
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Figure 18. 
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Torsion creep testing setup in a different location with 
furnace door open. 
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Figure 19. Troptometer bridge circuit. 

E. Dry cell, 1V 2 V. 

M. Null indicator, Honeywell 104 WIG. 
R D . Decade resistance box, to 11,110 ft. 
R F . Precision fixed resistor, 5,000 ft. 
R L . Left potentiometer, to 10,000 ft. 
R R . Right potentiometer, to 10,000 ft. 



corresponds to an angle of twist of 2.42 X 10~ 6 radians 
per inch. 

Since an exact solution is available for elastic torsion 
of circular tubes, a torsion test of a circular tube 
specimen was carried out at room temperature, within 
the elastic range, to verify the loading and twist 
measuring systems. The shear modulus derived from 
this test was 3.98 X 10 6 lb/in 2 which agrees with the 
nominal value of 4.0 X 10 6 lb/in 2 for 2024-T3 aluminum 
alloy. 

C.2. Procedure 

Twelve torsion creep tests were carried out at 400 
°F, three on each of the four tube shapes shown in 
figure 1. 

Alinement of each specimen was visually achieved 
by adjusting the positions of the bearings which support 
the torque shaft. (This procedure was verified with 
one 0.083-in wall square tube specimen using resist- 
ance strain gages. Thirteen gages were cemented 
onto the outer surface of the specimen, at various 
locations, in the longitudinal direction. Torques were 
applied, within the elastic range, at room temperature. 
The strain measurements did not reveal the presence 
of any significant bending or axial force in the speci- 
men.) The weight pan was supported on a platform 
such that the sprocket chain was unloaded, and 
weights of the desired amount were stacked upon it. 
Specimens were heated to the test temperature, 400 
°F, in approximately 1 hr, and were kept at this tem- 
perature for another hour to permit temperatures to 
equalize prior to the start of the creep test. The test 
was started by smoothly raising the weight pan from 
its platform with a turnbuckle (N, fig. 17) and then 
removing the platform. 

Twist and temperature measurements were made 
at closely spaced intervals following the start of a 
test and less frequently thereafter. Twist readings 
were found to be repeatable within 1 11 on the decade 
box (which is equivalent to an angle of twist of less 
than 3 X 10 -6 radians per inch) and temperatures were 
constant within 1 °F throughout the tests and uniform 
within 3 °F. 

'Tests were terminated after approximately 50 hr. 

Appendix D. Anisotropic Torsion Creep Theory 
for Noncircular Tubes 



length. A bridge circuit, figure 19, is used to measure 
the difference in resistance between the left and right 
potentiometers. The decade resistance box (L, fig. 
16 and R D , fig. 19) is periodically adjusted to give a 
null indication on the meter (M, figs. 16 and 19) and 
the box setting is then equal to the desired resistance 
difference. The sensitivity of the null indicator is 
equivalent to a fraction of an ohm and resistance 
measurements were made to the nearest ohm. With 
the sprocket ratios and potentiometers used, one ohm 



Effects of anisotropy in plastic deformation are 
commonly included by using a modified expression for 
the octahedral shear stress, such as 



_ 1 



t = 3 [ Wcr* ~ <ry? + «cr„ - cr z f + U<r* ~ <^) 2 + 6£ 4 r; 
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+ 6^ + 6&rf x ]'/ 2 (D.l) 

where x, y, and z are Cartesian coordinates which 
denote the principal axes of orthotropy and & are con- 
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stants. Coordinate x is parallel to the longitudinal axis 
of the tube. In this form r does not have physical sig- 
nificance as the shear stress on the octahedral plane 
unless all the £j are unity; it is simply a measure of the 
propensity for plastic flow (in this case creep). The 
octahedral shear creep strain is, as usual, 



■=f[(€*-€ y > 



+ (€y-€ g ) 2 + (€ g -€ x ) 2 
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Solving eqs (D.3) and (D.4) simultaneously and substi- 
tuting eqs (1) and (3) therein gives 



where 
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\/ 7~T"7 is the amsotropy parameter. 
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Performing a contour integration to evaluate y m and 
Using eqs (D.l) and (D.2) it may be seen that in uniaxial substituting into eq (2) gives, for a tube of uniform wall 
tension with a x = a and e x = e, and assuming Poisson's thickness 
ratio is one-half, 



7= VCl + 6, (randy=V2e. 
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c . .i i r i tu — ™^ ^, — », When Z=l, this agrees with the isotropic twist equa- 

Similarly, for pure shear with r xy — T and y X y— 7 .. ,-v' j u y lie •♦ *u \u 

7 tion, eq (7), and when Z = 1.15 it agrees with the 
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anisotropic version, eq (10). 
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(Paper 71C3-256) 
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